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A STUDY OF THE WAKE AT SUPERSONIC AIR FLOW 
OF AN ELLIPTICAL BODY - FLAT DISK COMBINATION 
USING HYDRAULIC ANALOGY TO COMPRESSIBLE AIR FLOW 


SUMMARY 

The problem investigated herein is the Mach number recovery 
in the wake of an elliptical cylinder in a supersonic flow whose 
free stream Mach number is approximately 2.00. Using the hydraulic 
analogy to simulate conditions in compressible gas flow, Mach 
number recovery is determined by measurement of the distances be- 
tween the detached shock waves and a flat disk model towed at 
various distances behind the elliptical cylinder model, 

Comparisons are made between the theoretical forms and loca- 
tions of detached shock waves ahead of blunt bodies in air and 
the observed forms and locations of detached hydraulic jumps ahead 
of blunt bodies in water, 

Condensed hydraulic analogy and detached shock wave theories 
are included in the appendices. 

As a result of this investigation, the following conclusions 
are made: 

a. A water channel provides a very inexpensive method 
for demonstrating two-dimensional wave characteristics in com- 
pressible gas flow. Within the limitation of the hydraulic anal- 
Ogy, quantitative data of the wave formation can be obtained by 
carefully controlling the water depth and velocity. 

b. The hydraulic analogy is valid for use in the study 
of leading edge shock waves caused by flow deflections approaching 


the shock wave detachment angle for a double-wedge, 


c. The detached shock wave theory very closely predicts 
the forms of detached shock waves ahead of blunt bodies, 

d. Mach number recovery in the wake of an elliptical 
cylinder at a free stream Mach number of approximately 2,00 is 
90 percent complete at a downstream distance of approximately 2.4 
model widths, Recovery is 95 percent complete at a downstream 
distance of approximately 3.4 model widths. 

e, At downstream distances equal to or greater than 
seven times the model width, the flow between the trailing edge 
shock waves of an elliptical cylinder is uniform, 

f, Until confirmed by experimental evidence obtained in 
a supersonic wind tunnel, the quantitative data in d above must be 


viewed only as a guide for use in planning wind tunnel experiments, 





A STUDY OF THE WAKE AT SUPERSONIC AIR FLOW 
OF AN ELLIPTICAL BODY - FLAT DISK COMBINATION 
USING HYDRAULIC ANALOGY TO COMPRESSIBLE AIR FLOW 


INTRODUCTION 

This investigation was conducted for the purpose of determin- 
ing the Mach number in the wake of an elliptical cylinder at down- 
stream distances up to ten times the width of the cylinder and the 
observation of wave formations caused by various models, A free 
stream Mach number of approximately 2.00 was maintained during 
all test runs. 

The hydraulic analogy to compressible gas flow has been 
demonstrated by Preiswerk(1)* and others. The simulation of shock 
waves in air by hydraulic jumps in water is the basis upon which 
this investigation was conducted. The form and the location of 
the detached shock waves produced by a flat disk in the wake of 
the elliptical cylinder provide the quantitative data for this 
report. 

This experinent was conducted at the University of Minnesota, 
Minneapolis, Minnesota, during the school year 1956-1957. 

It is desired here to express thanks to some of those whose 
assistance made this experinent possible; namely: to Professor 
John D, Akerman and Dr. Helmut G. Heinrich for their inspiration 


* Numbers in parenthesis refer to references listed at the end of 
this report, 








and many pertinent suggestions; to Mike Schonberg and Miles Mock 
for providing the necessary hardware; to Arvid Grikis for his 
assistance with the photographs and for reproducing the report; 

and to Stanley Markham and John Kallevig for the myriad jobs under- 
taken in the conduct of the investigation. Finally, to my wife, 
Marty, is extended my most sincere thanks for her uncomplaining 
endurance of three years of virtual solitude during the current 


"shore duty" tour. 
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TABLE OF SYMBOLS 
speed of sound in a gas 
cotangent of the Mach wave angle ( Yv2 — 1) 
specific heat at constant pressure 
specific heat at constant volume 


adiabatic gas constant, Cp/cy 


water depth 

angle between sonic line and normal to free-stream direction 
acceleration due to gravity 

enthalpy or novall heat content 

angle of streamline relative to x-axis 

wedge half-angle for which shock becomes detached 

Mach number, stream values unless otherwise indicated (V/a for 
gas; V/ Ved for water) 

Mach wave angle (sin71 =) 

pressure 

Mass density 

isentropic contraction ratio from free stream to sonic 
velocity 

absolute temperature of gas 

velocity 


components of velocity in x-direction and y-direction, 


respectively 


X3y,Z rectangular coordinate axes 








SUBSCRIPTS 
centroid of stream tube passing sonic line 
value at stagnation (V = 0) 
value in undisturbed stream 
sonic point of detached shock 
sonic point on body 
conditions along sonic line 
conditions in the wake of the elliptical cylinder model 


partial derivative with respect to x (for example 


bx. ~ 28, bsex is 2b) 


partial derivative with respect to y 














EOQULFMENT 

Three models were used during this investigation, a photo- 
graph of which is included as Fig. 1. They are: 

a. Double-wedge - constructed of pine lumber, 6.00 
inches in length, having a half-angle of 20 degrees, waterproofed 
with boiled linseed oil and varnished to provide a smooth surface. 
The height of this model is approximately lt inches. 

b. Elliptical cylinder - constructed of pine lumber, 
2.00 inches in length, 1.00 inch in width, approximately 12 
inches in height, treated in the same manner as the double-wedze. 

c. Flat disk - constructed of brass sheet stock 0.040 
inch thickness, 1.00 inch wide and approximately one inch in 
height. A thin brass rod is soldered to this model to facilitate 
its being mounted onto the model tow carriage, 

The models may be attached to tubular arms which extend out 
perpendicularly from model carriages mounted on an extruded 
aluminum monorail. This monorail is suspended approximately eight 
inches above the plate glass floor of the water basin used in this 
experiment. A photograph of the model set up is shown in Fig. 2. 

The model carriage towing system is composed of a power and 
velocity control unit for the tow, and a closed circuit phosphorus- 
bronze stranded wire towline. The power and control unit is a 
115 volt A. ©. variable speed, reversible motor producing up to 
five inch-pounds of torque with an electronic speed adjustment, 
This unit is shown in Fig. 3. 

The water basin used in this experiment is L-shaped, is four 


feet wide and has a maximum length of twelve feet. The plate 
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glass floor of the basin may be leveled by adjusting the height 
of the supporting legs. A complete description of the design of 
this basin is found in Ref. 2, 

In order that the model speed could be measured accurately, 
a timing mechanism was constructed, consisting of a thin, formed 
copper contact, mounted on the forward model carriage, which is 
held by spring action against two thin brass rails inserted in 
grooves in a wooden beam mounted alongside the model carriage 
monorail. One of the brass rails is 38.8 inches in length,is 
mounted opposite the test section and serves as the timer strip, 
while the other rail extends over the entire length of the wooden 
beam. In order that the copper contact brush may smoothly trav- 
erse the entire length of the water besin, plexiglass strips are 
placed in the grooves on either side of the timer brass strip. A 
photograph of the timer mechanism is shown in Fig. 4. Each brass 
rail is connected by wiring to a modified Brush Oscillograph 
Recorder control unit (Fig. 5) in such a way that when the copper 
contact brush closes the circuit between the two brass rails, a 
60-cycle signal is introduced to the oscillograph. When the model 
carriage has passed the test section and the copper contact brush 
has broken contact with the timer rail, this signal is abruptly 
interrupted. Paper tape is fed from the oscillograph recorder at 
the rate of 125 millimeters per second. Since the model speed is 
low (of the order of 1.6 feet per second), this timing device is 
quite accurate, 

Since the speed of propagation of gravity waves in liquids is 


dependent upon the depth of the liquid (see Appendix A), it is 












necessary to measure the depth of the liquid accurately. For this 
purpose a device was constructed, eo of a micrometer 
mounted in a framework of aluminum sheets and wooden blocks, con- 
nected by wiring to an ohmmeter containing a small battery. A 
photograph of this device is included as Fig. 6. Tap water used 
in the water busin contains a sufficient amount of ionized minerals 
to be a good electrical conductor, One electrical lead is con- 
nected to the micrometer head and the other lead is placed in 

the water. When the micrometer shaft contacts the surface of the 
water, the ohmmeter needle is deflected. With this device, water 
depth can be measured with O0,0Q0l~inch accuracy. Caution must be 
exercised, however, in insuring that the micrometer shaft is 
perfectly dry and free of foreign matter in order that premature 
contact with the surface of the water is avoided, 

With the use of the timer mechanism and the water depth 
measuring device, it is possible to compute the test section Mach 
number with accuracy to the second decimal place, 

Photographs were made with a studio type camera mounted 
vertically on a tripod supported above the test section by a wooden 
frame, the camera lens being approximately four feet above the 
water line, Exposure times for all test run photographs were of 
1/50 second duration, Lighting was provided by an electronic 
flash unit, the duration of the flash being 1/500 second, For the 
test runs, aperture openings of f16 and f22 were used. Eastman 
Super~Panchro, Type B, 4" X 5" sheet film was used to make all 
photographs, 


Several photographic lighting arrangements were tried. It 





wis found that for the observation of the detached shock waves aheed 
of the models, the most desirable lighting arrangement is obtained 
when the electronic flash is placed slightly above the camera, é&imed 
at the surface of the water slightly upstream of the models, The 
light reflected from the disturbed surface of the water provides a 
very good high-contrast photograph. A black velvet cloth was sup- 
ported beneath the glass floor of the water basin to provide a dark 
background for the photographs. A view of the entire experimental 


set up is included as Fig. 7. 


PROCEDURE 

In what follows, a general knowledge of the theory of the hy- 
draulic analogy to two-dimensional compressible air flow(1) is 
assumed, A summary of this theory is included as Appendix A, 
Appendix B contains a summary of an approximate method for predicting 
form and location of detached shock waves ahead of plane bodies(3) , 

Since Laitone 4) has analyzed the theory governing the propa- 
gation of surface water waves to show that an analogy with the tio- 
dimensional flow of a perfect gas exists only if the water dept! is 
approximately one-fourth inch, water depths in this experiment were 
maintained as near to 0.250 inch as practicable, 

It is also known that the hydraulic analogy to compressible air 
flow becomes less valid at high Mach numbers. It was, therefore, 
decided to use a rather low supersonic Mach number for these test 
runs, Mach amniger 2.00 was arbitrarily chosen, 

The following series of test runs were made: 


a. Double~wedge alone 





b. Elliptical cylinder alone 

c. Flat disk alone 

d. Flat disk placed directly in the wake of the elliptical 
cylinder at distances varying in one-inch increments from one inch to 
ten inches behind the trailing edge of the elliptical cylinder, 

The bottom surface of each model was in sliding contact with the 
water basin floor at all times during each run. 

The general procedure in making all the test runs was as 
follows: 

a. Check and align model(s) to insure that the model(s) 
had no angle of attack, 

b. Adjust depth of water to be as near 0,250 inch as 
practicable. 

c. Using the timing device, make dummy runs, adjusting 
the motor control unit to provide model speed as near desired sim- 
ulated Mach number 2.00 as practicable. 

d, Check camera focus, using a printed cerd floating on 
the surface of the water as a focus target, 

e. Again check and record the depth of the water. 

f. Make test run, taking the photograrh as the model (s) 
pass through the test section, | 

g. Record photographic data and test run Mach number, 

Due to low relative humidity in the room containing the water 
basin, measurement of water depth had to be made very shortly 
before making each test run in order that the correct Mach number 


could be computed, 








10 
RESULTS AND DISCUSSION 

The results of this investigation are contained in Table I and 
in Figs. 8 through 26. Fig. 8 is an overlay, depicting the theoret- 
ical shock wave pattern associated with a 20-degree half-angle double- 
wedge model in air. In Fig. 9, the actual shock wave produced by 
the double-wedge model in water is shown. The theoretical detached 
shock wave pattern ahead of the elliptical cylinder are shown in an 
overlay in Fig. 10. Fig. 11 shows the actual shock wave pattern 
associated with the elliptical cylinder. Fig. 12 is another over- 
lay, showing the form and location of the theoretical detached shock 
wave ahead of the flat disk. Fig. 13 shows the actual shock wave 
associated with the flat disk. Based on the results of comparisons 
made in Figs, 12 and 13, a modified detached shock wave location 
theory for the flat disk model is plotted in Fig. 14. 

In Figs. 15 through 20 and Figs. 22 through 25, shock wave 
patterns associated with the combination of the ellipticel cylinder 
and the flat disk models are shown, Model separation distances in 
one~inch increments from one to ten inches downstream of the trailing 
edge of the elliptical cylinder model are shown. The theoretical 
detached shock wave ahead of the flat disk model is shown as an 
Svemlay an Fig. 21. 

Mach number recovery in the wake of the elliptical cylinder is 
tabulated in Table I and plotted versus downemmermete canes an 
eile ee 

It is desired to establish the validity of comparing the form 
of hydraulic jumps in water caused by large flow deflection with 


shock waves produced under similar conditions in air. A double- 





tal 


wedge model, upon which is impinged an attached leading edge shock 
wave, 1S ideal for making such comparisons. Theoretical shock weve 
engles for double-wedges in air flow are ecsily computed from Ref. 6. 
The wedge half-angle reaouired for shock wave detachment in air at 
Moo = 2.00 is approximately 23 degrees. Therefore, a model whose half-~- 
angle is 20 degrees provides insurance that the shock wave will be 
attached even though the flow has been deflected through a signifi- 
cantly large angle. It is pertinent here to note that the leading 
edze shock weve remains straight only until it enccunters the Mach 
w:.ves caused by flow expansion around the model at its maximum width. 
For the model used in this exreriment, this interaction takes place 
(at My = 2.00) approximately three inches in the iomcigae direction 
from the leading edge of the model. 

A comparison of the theoreticel and measured shock wave angles 
associated with the double~wedge model are shovm in the unnumbered 


table below, 


Theoretical shock wave Measured shock 

angle (¥=1.4) wave angle 
Leading edge shock 53,2 544° 
Trailing edge shock 20, 17,5° 


it may be seen that the observed leading edge shock wave angle is in 
Quite close agreement with the compressible gas shock wave theory. 
The discrepancy between theoreticel and measured trailing edge shock 
wave angles is significant (approximately 15 percent). This discrep- 
ancy may be partially explained by the fact that the hydraulic 
snalogy becomes less valid for flow deflections at higher Mach 


niibers. The theoretical Mach number on the rear fece of the dontle- 
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wedge model is 2.7k. 

From the foregoing, it may be deduced that the hydraulic analogy 
is applicable, with certain reservations noted below, for the study 
of shock wave formations caused by large flow deflections, 

In Figs. 10 and 11, comparison is made between the forms and 
locations of theoretical and observed detached shock waves ahead 
of the elliptical cylinder model. The dimension yop (see Appendix 
B) for an elliptical body is determined by the shock weve detach- 
nent angle, which is dependent upon the adiabatic constant v. 
Therefore, theoretical shock wave forms and locetions for both air 
(¥ = 1.40) and the analogous gas (¥8= 2.00) are shown. It is seen 
immediately that the locations of the theoretical shock waves are 
at considerable variance with the location of the observed shock weve. 
However, the forms of the theoretical shock wave for the analogous 
gos (¥= 2.00) and the observed shock wave are almost identical. It 
may be concluded that the theory developed by Moeckel (3) very closely 
predicts the form of a detached shock wave ahead of an elliptical 
cylinder model in water. However, as is also shown in Ref. 3, ex- 
perimental data indicete that the theory places the location of the 
detached shock wave too close ahead of a blunt-nosed body in air, 

Similar conclusions to those mentioned above may be drawn con-~ 
cerning the detached shock wave ahead of the flat disk model. In 
Fig. 12, a theoretical detached shock weve for air (¥8=1.4) is 
shown, The form of the theoretical shock wave in air was found to 
be very nearly the same as the observed form of the hydraulic jump 
in water. In view of the patterns shown in Fig. 10, it appears 


reasonable to assume that the form of the theoretical shock wave in 








ae 


the analogous gas would be virtually identical to that of the observed 
pattern in water. 

In order that the observed position of the detached shock wave 
ahead of the flat disk may be used to determine the Mach number 
ahead of the model it is necessary to modify the theory in the 
Vicinity of Mg= 2.00. The distance L from the model to its de- 
tached shock wave (Fig. 13) at Mg= 2.00 is measured to be 0.85 inch. 
Since ysp for the flat disk model is 0.50 inch, the dimension 
L/yggp is 1.70. The dimension L/ysp for Mg -= 2.00 in accordance 
with Ref. 3 is 2.12. The ordinates of the theoretical curve of 
L/ygp versus M, from Ref. 3 are reduced by the ratio of 1.70 
to 2.12 in Fig. 14. This reduction appears to be justified, in view 
of the fact that all measured data are in the vicinity of Mg= 2.00. 
Further, experimental data recorded in Ref. 3 indicate that it is 
to be expected that detached shock waves ahead of very blunt bodies 
will be farther removed than is predicted by the theory. 

The location of the detached shock wave ahead of the flat disk 
in the wake of the elliptical cylinder is not perfectly clear in all 
photographs. Therefore, the measurement of the dimension L becomes 
largely a matter of judgment. In Figs. 17 through 20, the assumed 
locations of the foremost points of the detached shock waves are 
indicated. 

The distance L from the flat disk model to its associated 
detached shock wave for each test run is tabulated in Table I. The 
value of L/ysp is then used with Fig. 14 to determine the corre- 
sponding Mach number. To make the Mach number recovery presentation 


nondimensional, each measured Mach number is divided by Mg. A plot 
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Pa 


of My /Mop versus downstream distance from the elliptical cylinder 
model is made in Fig. 26, using the shock wave positions as positions 
of measured wake Mach numbers. Due to the rather large distances 
between the shock waves and the model, it is believed that the po- 
sitions of the shock waves themselves more nearly define the posi- 
tions of measured Mach numbers than do the positions of the model. 

It is noted here that for model separations of one and two inches 
Figs. 15 and 16), no characteristic detached shock wave appears ahead 
of the flat disk model. It must therefore be concluded that the 
flow immediately behind the trailing edge of the elliptical cylinder 
model is subsonic. This, of course, is to be expected in viscous 
flow. 

It is seen in Fig. 26 that the wake Mach number recovery in the 
analogous gas is 90 percent complete at a downstream distance of ap- 
proximately 2.4 model widths and is 95 percent complete at a down- 
stream distance of approximately 3.4 model widths. Beyond that 
point, the wake Mach number recovery is quite gradual, approaching 
100 percent recovery asymptotically. 

The first downstream flat disk model position for which a com- 
plete, clearly defined detached shock wave is visible between the 
trailing edge shock waves of the elliptical model is shown in Fig. 
22. Since it is desired to compare the form of the detached shock 
wave in this area with the theoretical shock wave, an overlay is in- 
cluded as Fig. 21. In Fig. 21, the theoretical detached shock wave 
for the flat disk in a free air stream corresponding to the computed 
wake Mach number is shown. This theoretical shock wave is corrected 


for location in accordance with the modified L/ysgp curve shown in 
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Fig. 14. It can be seen that, except for the deviation caused by 
interaction between the detached shock wave and the trailing edge 
shock waves of the elliptical cylinder, the observed form of the 
hydraulic jump agrees fairly well with the theory in a free air 
stream. In view of the comparison made between the shock waves in 
air (8 = 1.4) and the analogous gas (¥= 2.0) in Fig. 10, it is 
suspected that the ieeeerieed detached shock wave for the analogous 
gas would be very nearly identical to the hydraulic jump observed 
in Fig. 22. A similar comparison also was made with Fig. 25, the 
results being the same. It may be concluded, therefore, that at 
downstream distances equal to or greater than seven times the model 
width, the flow between the trailing edge shock waves of the 
elliptical cylinder is uniform, 

A search of the available literature fails to reveal a previous 
investigation of this nature in air. Therefore, the data provided 
by this experiment are not corroborated. In view of the foregoing, 
it must be concluded that until confirmed by investigators using two- 
dimensional models in a supersonic wind tunnel with a Reynolds number 
of approximately 22,500 based on the length of the elliptical cylinder 
model, the quantitative data shown herein must be viewed only as a 


convenient guide for use in planning wind tunnel exreriments, 





' 
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CONCLUSIONS 
From the results of this investigation the following conclusions 
may be made, 

a. <A water channel provides a very inexpensive method for 
qualitatively demonstrating two-dimensional compressible gas flow 
phenomena. Within the limitations of the hydraulic analogy, quan- 
titative data can be obtained by carefully controlling the water 
depth and velocity. 

b. The hydraulic analogy is valid for use in the study of 
leading edge shock waves caused by large flow deflections, 

c. The theory develored by Moecke1 (3) very accurately pre- 
dicts the forms of detached shock waves ahead of blunt bodies. How- 
ever, as is also shown in Ref. 3, this investigation indicates that 
the theory places the locations of the detached shock waves too 
close ahead of blunt bodies, 

d, Mach number recovery in the wake of an elliptical cyl- 
inder at My= 2.00 is 90 percent complete at a downstream distance 
of approximately 2.4 model widths. Recovery is 95 percent complete 
at a downstream distance of approximately 3.4 model widths. 

e. At downstream distances equal to or greater than seven 
times the model width, the flow between the trailing edge shock wuves 
Seecnel elliptical cylinder is uniform, 

f. Until confirmed by experimental evidence obtained in a 
supersonic wind tunnel, the quantitative data shown herein must be 
viewed only as a convenient guide for use in the planning of wind 


tunnel experiments. 
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MACH NUMBER IECOVERY IN THE WAsE OF AN &LLIPTICAL CYLIUDER 
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APPENDIX A 
SUMMARY OF HYDRAULIC ANALOGY THEORY 

It has been known for some time that an analogy exists between 
the flow of a liquid with a free surface and the flow of a compress- 
ible gas. Early mathematical treatments of the analogy were made 
by Jouget and Riabouchinsky. In 1938, Preiswerk(1) made the first 
really comprehensive theoretical and experimental investigation of 
this phenomenon. In more recent years, a great deal of interest has 
been shown by numerous groups in the hydraulic analogy. Especially 
notable among these groups has been the National Advisory Committee 
bor Aeronautics(5), 

The following is a summary of the theory and mathematical develop- 
ment of the hydraulic analogy. 

Two assumptions are made in the mathematical development: 

a. The flow is irrotational. 

b. The vertical accelerations at the free surface of the 
liquid are negligible compared with the acceleration of gravity. 
Therefore, the pressure in the liquid at any point depends oniy on 
the height of the free surface above that point. 

Consider a water channel consisting of a large tank in which the 
depth of the water, d,, is constant, and an open duct leading from 
the tank in which the depth of the water, d, varies with distance 
from the tank. The floor of the open duct is an extension of the 
floor of the tank, and is perpendicular to the gravity field of the 
Rarth, lf the vertical coordinate, z, is measured from the floor of 


the channel, the Bernoulli equation is: 


p+ Ave + Pez = Po + PB2 
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or 

V2 = 2¢(z, - z) + 24Po - P) (A-1) 
Utilizing assumption b above: 

Po = gldg - &) (A-2a) 
and 

p = g(d - z) (A-2b ) 
Direct substitution into equation (A-1) gives 

V“ = 2e(dg - d) = gad (A-3) 
Since equation (A-3) does not contain z, the velocity at any point 
in the open duct is constant over the entire depth and is given only 
by the difference in depth Ad between the total head and the free 
level. Ad, at the most, is equal to d,. Therefore, the maximum 
attainable velocity is 

Vmax = ¥2edo (A-L) 
The corresponding equations for the velocity of a gas are 

V* = 2e(ho -h) = 2gep(To - T) (A-5a) 
and 


Vmax = Y28No = ¥2gepl (A=5b) 


Therefore, if the ratio V/Vmax for gas is equated to V/Vyax for 


water, 
Gn = 6l 20h Sri Ws <a at 
ae ie To 
or 
a =T 
do To (A-6) 


The continuity equation expresses the fact that, for a prizam of 
Space in an incompressible fluid, outgoing volume equals the in- 


coming volume . The equation of continuity for water is 





LO 


Olud) , _Olvd) _ 9 
Ox Oy (A-7) 


The continuity equation for two-dimensional compressible gas flow is 


SON) perol2))) 
Oo x Oy 


(A-8 ) 
It may be seen that equations (A-7) and (A-@) are of the same form. 
From these equations we may derive a further condition for the analogy, 


that the density of gas flow corresponds to the water depth d. Thus, 


ol 
aS (A-9) 


Previously, however, it was seen that water depth d was the 
analogous magnitude for the gas temperature T (equation (A-4), 
This is possible only if a very definite assumption is made regards 
the nature of the comparison ges. For adiabatic isentropic flow in 
a gas 


! 


ao 


and since, in the analogy, from equations (A-6) and (A-9) 


a 
~ % Te 


! 
Dee 
To To 


which is satisfied only if ¥ = 2.0. 


then 


From the gas relation 


it may be seen, since Y= 2.0, that 





Al 


i ' (2) (A-10) 


The velocity potential for water is given by the equation 


ter (1 - 22) + by (2 - 8) ~ 2bgy bY 2 (aan) 


x 
gd gd 


and the corresponding equation for a gas is 


sex (2 - ts ty [2 - ) - 2bxy hy = 0 (A-12) 





Therefore, for identical expressions, 
gd = af (A-13) 

From equations (A-11) and (A~-12), the velocity Vgd in the 
liquid flow is seen to correspond to the velocity of sound in gas 
flow. It has been shown (5) that Ygd is the velocity of propagation 
of surface waves, the wave lengths of which are large in comparison 
with the water depth. The ratio of V/ gd in the liquid flow cor- 
responds to the Mach number V/a in the gas flow. 

If the velocity of the liquid flow is less than Ygd (M<1), 
the water is said to be "streaming". If the velocity of the liquid 
flow is greater than Ved (M>1), the water is said to be "shooting". 

In shooting water under certain conditions, the velocity of the 
flow may strongly decrease for short distances and the depth may 
increase, An unsteady motion of this type is called a hydraulic 
jump. Hydraulic jumps of small intensity are propagated with the 
velocity Ved, corresponding to Mach waves in compressible gas flow. 


The analogy may be summarized in tabular form as follows: 





Significant quantities and 
characteristics of two-dimensional 
compressible gas flow, %= 2.0 


Temperature ratio, T/T, 
Density ratio, Gs 


Pressure ratio, p/P 


Velocity of sound, a = -e 
fa 


Mach number, V/a 


Subsonic flow 
Supersonic flow 


Shock wave 
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Corresponding values in 
analogoug liquid flow 


Water-depth ratio, d/dy 
Water-depth ratio, a/do 


Square of water-depth 
ratio, (d/do)2 ‘ 


Wave velocity, Ved 


Mach number, V/ gd 
Streaming water 
Shooting water 


Hydraulic jump 
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APPENDIX B 
FORM AND LOCATION OF DETACHED SHOCK 
WAVES AHEAD OF PLANE BODIES 

The following is a condensation of pertinent portions of Moeckel's 
approximate method for predicting form and location of detached shock 
waves ahead of plane bodies (3). Computations used in predicting the 
form and location of detached shock waves ahead of models used in 
this experiment are included in Appendix C. 

The representation of flow with detached shock waves used in 
this analysis is shown in Fig. B-l. The region of subsonic flow is 
bounded by the portion of the detached wave below its sonic point, 

S, by the portion of the body contour upstream of its sonic point 

SB, and by the sonic curve between S and SB. If the form of the 
detached wave is known, the location of S on the shock can im- 
mediately be determined, because for a given Mg the angle Ws for 
which sonic velocity exists behind the shock is known from shock 
Preorao) The flow direction As at this point is then also known. 
The sonic curve is shown in Fig. B-l as a straight line, This ap- 
proximation is used throughout the analysis, although more exact com- 
putations, as well as experimental results, indicate that the form 
may depend considerably on the shape of the nose or leading edge. 

To the degree of approximation of the present method, however, such 
variations appear to be unimportant. With the simplified picture 
shown in Fig. B-l, approximate expressions can be derived for the 
shock location relative to the body sonic point. 

Evidence is available to show that, for bodies with sharp or 


well-defined shoulders, the sonic point is located at the shoulder. 





ae | 


Aly 


For more gradually curved bodies, such as ogives or ellipses, the 
location of the sonic point may be estimated to be at the point where 
the contour of the body is inclined at the wedge angle corresponding 
to shock detachment. 

The most typical characteristics of detached waves are that: 
(a) they are normal to the free stream at their foremost point; and 
(b) they are asymptotic to the free-stream Mach lines at large 
distances from their foremost point. A simple curve that has these 


characteristics is an hyperbola represented by 


By =\x* - x0* (B-1) 

where 6 is the cotangent of the Mach angle and xo is the distance 
from the vertex of the wave to the intersection of its asymptotes. 
For the purpose of the analysis, the hypothesis is made that to the 
degree of approximation required, all detached waves in the region 
between the axis and the sonic point S may be represented by equation 
(B-1). 

With this form of detached wave, the angle between the stream 
direction and the tangent to the shock at any point is obtained from 


(Fig. B-1) 

dy = tan Y = x = 1x0? + Py? 

dx A x2 = x02 8 -y (B-2) 
The location of S is then 


xQ cot Ps 


a 6 |p? - cot UP (B-3 ) 


cog Oe 
Ve cot S (Bet, ) 


If the y-coordinate of the body sonic point is used as the reference 


dimension, then the form of the shock wave is given by 
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Vv a oe x \2 _ XO Z 
Ysp_ so Bs aa a (B-5) 


where, from equation (B-3), 


ne 2 7 
YSB <= | 2 ee - 2 (B-6) 


From the equation (B-4) the dimensionless location of the shock sonic 


point is 
X0 
28 . Weise 
YSB V6 eee (B-7) 


The distance from the foremost point of the shock to the x-coordinate 


of the body sonic point is 


+ =x9B - % 
YSB YSB YSB (B-8) 


where, from Fig. B-l, 


*SB _ 4S (3 


-~ 1} t Be 
YSB YSB SB eon — 


If equations (B-6), (B-7) and (B-9) are combined to eliminate all 


unknown coordinates except yg and L, equation (B-8) becomes 


au 2 ee (C+tan eon y {B-10) 
Ysa Ysp 


where 
fp (Bran Y = I a * tan“ Y - 1) (B-11) 
Inasmuch as p and Un are known for any given free-stream Mach 


number, only quantities ys/ysp and " remain to be determined to 


predict the relation between the sonic point on a body and the 
location of its detached wave, 
The continuity equation is applied to the fluid that passes the 


sonic line to provide 





1,6 


om = E = (Po )eo 
Po/e 


go cos n | . 


(B-12) 
where o is the contraction ratio required to decelerate the free 
stream to sonic velocity isentropically. 

The appropriate value of y to be used remains to be established, 
On the basis of the analogy between one-dimensional channel flow and 
flow with detached shock waves, the sonic line is assumed to be normal 
to the average flow direction in its vicinity. At 5S the incli- 
nation of the flow is known to be A 5; whereas at SB the inclination 
is assumed to be Nae If the arithmetic mean of the inclinations 


at the two extremities is used, the expression for n becomes: 


n =% ( Na + As) (B13) 


Because Ng differs only slightly from Ag, the inclination of the 


sonic line for plane flow is assumed to be simply = As, 


Values of Ys/Ysp obtained from equation (B-12) are shown as 


functions of My in the neighborhood of Mg = 2.00 in Fig. B-2. The 
variation of Xo/yon with My, as computed from equation (B-6), is 


shown in Fig. B=3. Values of L/Yop obtained from equation (B-10) 


are shown in Fig. Bed, 
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Fig. B-3 
Parameter used for determination 
of shock wave form and location 
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Fig. Be4 
Theoretical shock wave location 
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APPENDIX C 
SAMPLE CALCULATIONS 
1. Predicted form and location of detached shock waves ahead of 
models. 

It is seen in Appendix B that the reference dimension upon which 
to base the prediction of the form and location of the detached shock 
waves ahead of a model is Ysp. This dimension for an elliptical 
model, is dependent upon the shock detachment angle Xa for a wedge 


in two-dimensional flow. From Ref. 6 


» E I (me 2 = 1)1-? 
2 (Meo “ = 1)"*" 
te a Moo * (C-1) 


Substituting the analogous gas constant, Y= 2,00, and the test Mach 
number, M = 2,00, into equation (C-1l) gives 
Ag = 0.333 radians = 19.1 degrees 


The equation for the cross section of the elliptical model is 


ss yo =] 

DaaK(z)° (C~2) 
By setting the derivative of equation (C-2) equal to tan Ag, the 
x-coordinate of the body sonic point, measured from mid-chord, is 

Me= 2055075 ach 
Substituting this value into equation (C-2) gives 

eg = Opal inch (C~3 ) 
Measured from the leading edge of the model, the x-coordinate of the 
body sonic point is 

x = 1 = 0.5675 = 0.4325 inch 

For the flat disk model, the body sonic point is at the shoulder 


of the model, Thus, for the flat disk model 





Yop = 0-500 inch 
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(C-i, ) 


From Figs. B~2, B-3 and B-4, the following non-dimensional quan- 


tivues,sior ellie = 2.00 vars. found: 


ae 
YSB 

STE a 
YSB 

a0 6-96.65 
YSB 


From the foregoing, the following table may be computed: 


Quantity 
I 


YS 


XO 


XsB = x9 + L 


Elliptical Model 


1.480 inches 
7,670 inches 


8.542 inches 


Flat Disk liodel 


1.06 inches 
1.80 inches 


9,352 inches 


10,385 inches 


By substituting into equation (B-1), the theoretical shock wave 


coordinates for the two models are computed and plotted in Figs, C-l 


and G=2, 
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Theoretical shock wave for flat disk 
model in air (¥=1.4) at M,= 2.00 
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A study of the wake at super- 
sonic air flow of an elliptical 
body - flat disk combination 
using hydraulic analogy to com- 
pressible air flow. 





